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^ ■ Abstract 

O ■ 

i ■ Making use of the duality transformation, we derive in the Londons' limit of the Abelian 
Higgs Model string representation for the 't Hooft loop average defined on the string world-sheet, 
which yields the values of two coefficient functions parametrizing the bilocal correlator of the 
dual field strength tensors. The asymptotic behaviours of these functions agree with the ones 
Q H " obtained within the Method of Vacuum Correlators in QCD in the lowest order of perturbation 
theory. We demonstrate that the bilocal approximation to the Method of Vacuum Correlators is 
an exact result in the Londons' limit, i.e. all the higher cumulants in this limit vanish. We also 
show that at large distances, apart from the integration over metrics, the obtained string effective 
theory (which in this case reduces to the nonlinear massive axionic sigma model) coincides with 
the low-energy limit of the dual version of 4D compact QED, the so-called Universal Confining 
String Theory. We derive string tension of the Nambu-Goto term and the coupling constant of the 
rigidity term for the obtained string effective theory and demonstrate that the latter one is always 
negative, which means the stability of strings, while the positiveness of the former is confirmed 
by the present lattice data. These data enable us to find the Higgs boson charge and the vacuum 
expectation value of the Higgs field, which model QCD best of all. We also study dynamics of 
the weight factor of the obtained string representation for the 't Hooft average in the loop space. 
In conclusion, we obtain string representation for the partition function of the correlators of an 
arbitrary number of Higgs currents, by virtue of which we rederive the structure of the bilocal 
correlator of the dual field strength tensors, which yields the surface term in the string effective 
action. 
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1. Introduction 

Abelian Higgs Model (AHM) in the Londons' limit is known to display a lot of properties of the 
confining phase of gluodynamics, which makes it reasonable to consider Abrikosov-Nielsen-Olesen 
strings in this model as the objects which could tell us something about the behaviour of the 
gluodynamics string. To this end, in Ref. 1 the AHM partition function has been rewritten via 
string variables by making use of the duality transformation 2 and accounting for the Jacobian, 
which emerges when one passes from the collective field variables to the string ones. This Jacobian 
has been shown to yield the Polchinski-Strominger term 3 in the string action with the coupling 
constant ensuring exact cancellation of the conformal anomaly in 4D. 

An alternative approach to the string representation of gauge theories based on the Method 
of Vacuum Correlators (MVC) 4 ' 5 has been proposed in Ref. 6, where the QCD string effective 
action has been derived by performing the expansion of the Wilson loop average. Both in Refs. 1 
and 6, the sign of the rigidity term has been obtained to be negative, which means the stability of 
strings in the models under consideration 7 ' 8 . However, the two main disadvantages of the method 
proposed in Ref. 6 w.r.t. the one of Ref. 1 were the lack of the mechanism of cancellation of 
the conformal anomaly and the integration over metrics in the final expression for the Wilson 
average. The most natural way in the solution of these problems might lie in considering the 
strong background fields in QCD, which ensure confinement, as a collective variable leading to 
the summation over surfaces. 

In this Letter, we shall combine the approaches of Refs. 1 and 6 by considering the 't Hooft loop 
average in the AHM without monopoles and find the string representation for it. By making use of 
this representation, we shall demonstrate that the bilocal approximation to the MVC is an exact 
result in the Londons' limit, and obtain the values of two coefficient functions which parametrize 
the bilocal correlator of the dual field strength tensors. Studying the asymptotic behaviours of 
these functions at small and large distances and comparing them with the corresponding be- 
haviours known from the bilocal approximation to MVC in gluodynamics, we shall conclude that 
the Londons' limit reproduces well the long distance regime of gluodynamics, whereas the short 
distance one is reproduced only in the lowest order of perturbation theory. These points will be 
the topic of the next Section. 

In Section 3, we shall present the values of the string tension of the Nambu-Goto term and the 
rigidity term coupling constant in the obtained string effective theory. The latter one occurs to be 
always negative, which, as it was mentioned above, means the stability of strings. Demanding the 
positiveness of the string tension, we arrive at some constraint to the Higgs boson charge, which 
is satisfied according to the present lattice data. Based on these data, we also find the value of 
the v.e.v. of the Higgs field, which corresponds to gluodynamics best of all. 

In Section 4, we derive the equation of motion for the weight factor in the string representation 
of the 't Hooft loop average in the loop space. Both this equation and the string representation 
obtained in Section 2 demonstrate that, apart from the integration over metrics, the dual version 
of AHM in the Londons' limit coincides at large distances with the large distance regime of the 
dual version of the 4D compact QED, the so-called Universal Confining String Theory (UCST) 9,10 , 
whose loop space dynamics was investigated in Ref. 11. 

In Section 5, we obtain string representation for the partition function of the correlators of 
an arbitrary number of the Higgs currents, with the help of which we then get the correlator of 
two such currents. Due to the semiclassical connection between this correlator and the structure 
in the bilocal correlator of the dual field strength tensors which yields surface term in the string 
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effective action 5 , we rederive the coefficient function standing at this structure. 

In Conclusion, we shall focus ourselves at some physical consequences of the obtained results 
and discuss possible further developments of the suggested approach. 

2. String Representation for the 't Hooft Loop Average and the 
Exact Result for the Bilocal Correlator of the Dual Field Strength 
Tensors in the Londons' Limit 

In this Section, we shall consider string representation for the 't Hooft loop average, 
J Vx^) {J-'n(S)), which is defined on the closed string world-sheet S (we consider AHM without 
monopoles, where the only surfaces are the Abrikosov-Nielsen-Olesen strings' world-sheets). In 
what follows, we shall be interested in the weight factor of this representation, 




where $(x) = \$(x)\ e ie{x \ 6 = sing - + reg - , T^(x) = J da^(x(^))8(x-x^)) is the vorticity tensor 

s 

current, which is a functional of S, and = — iNeA^, so that the Higgs boson carries electric 
charge Ne. Notice, that, as it was already mentioned in the Introduction, in what follows, we 
shall, similarly to Ref. 8, for simplicity disregard the Jacobian arising when one passes from the 
integration over # smg - to the integration over £ M (£)- It is a l so worth mentioning, that from now on 
we shall make use of notations of Ref. 1, so that due to the lack of imaginary unit in front of the 
last term on the R.H.S. of Eq. (1), the sign of the second term of the cumulant expansion of Eq. 
(1) will be opposite to the standard one. 

In the Londons' limit, Eq. (1) takes the form 



(F N (S)) = J VA^VO^ exp^- J dx 
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(2) 



sing. 



Notice, that the S-dependence is present in Eq. (2) both in the last term and in the term <9 M # S 
due to the equation 2 e^ vp(T dpd a 6 snig \x) = 2tiT^ u {x). That is why, the naive application of the 
Gauss theorem (where only the S'-dependence in the last term in Eq. (2) is accounted for) is not 
valid, and therefore the S'-dependence in Eq. (2) could not be reduced to the volume one. 

Performing the duality transformation of the R.H.S. of Eq. (2) along the lines described in 
Refs. 1,2, and 8, we arrive at the following formula 



J VA^Vh^exp 



dx 



12r]' 



'iNe 7T 

— h^ + —T ap 



(3) 



where H pu \ = d^h v \ + d\h pu + d u h x ^ is a strength tensor of an antisymmetric field h^ v . From now 
on in this Section, we shall omit the preexponential factor, since our main result, the coefficient 
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functions parametrizing the bilocal correlator of the dual field strength tensors (see Eqs. (11) and 
(12) below) will be obtained upon comparison of the two exponents, the result of the Gaussian 
integrations and the one of the cumulant expansion of Eq. (2). 

One's first instinct is to fix some, for concreteness Feynman, gauge by adding to the Lagrangian 
standing in the exponent on the R.H.S. of Eq. (3) the gauge fixing term ^ (d^A^) 2 and carry out 
the Gaussian integration over A^, which then can be easily performed and yields 



J P^exp j-/ 



dx 



iNe 



-h. 



TV 



a/3 



exp I — - / dxdy 

\n 2 J \x- 



x-y\ 



4 ' Ne 

{x-y)^{x-y) u 
(x - y) 2 



- 1 )-\(W 2 



-5 



h^(x)h Xv (y) + 2niT Xtl (x)h Xv (y) 



(4) 



However we immediately see that Eqs. (3) and (4) provide us with such a representation of 
(Tn{S)), which is difficult to proceed with. 

It is much more convenient to rewrite exp dxF 2 ^ as 



/' 



VG^ exp 



dx 



after which A^- and G^-integrations yield 

(Fn(S)) = J VB^Vh^ exp jy dx 

2m 



(5) 



i.e. we have replaced the A^-integration by the one. Next, performing the hypergauge trans- 
formation with the function j^B^, which is equivalent to fixing the gauge by the condition B^ = 
(see the first Ref. in 8 ), we get from Eq. (5) 



{Fn{S)) = J Vh^exp jy dx 
which apart from the term 



u 2 



^) h 2 

2 J 



2 

- I TT- ) T 2 + 3mh uv T u 



(-)' 

\NeJ 



- jlV 



(6) 



exp 



(iVe) 



y dxT 2 u \ = exp | (j^j J da^ u (x) J do^x'^x - x' 



(7) 



is the partition function of the massive Kalb-Ramond field interacting with the string (nonlinear 
massive axionic sigma model 12 ). Notice, that this is also a low-energy expression for the partition 
function of the 4D UCST, studied in Ref. 10, which is in fact the dual version of the 4D compact 
QED. Therefore one may conclude that at sufficiently large distances, when the term (7) (or, better 
to say, its regularized version, which will be presented immediately below) can be neglected, the 
Londons' limit of AHM coincides with the low-energy limit of the 4D compact QED. 
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It is also worth mentioning, that the reduction of AHM in the Londons' limit to the nonlinear 
axionic sigma model enables one to apply methods developed for such models, for example, for 
the phenomeno logical description of fluctuations of the string world-sheet 13 . 

Regularizing the term (7) and performing the Gaussian integration over in Eq. (6) we 
arrive at the following representation for the weight factor (2) 



(Fn{S)) = exp 



(f)7 



)/ 



da^ v {x' 



-(x-x'f 



(8) 



-97r 2 y daxv(x) J da w {x')D Xv ^ p [x 
s s 

Here m = Net] is the mass of the Kalb-Ramond field (equal to the mass of the gauge boson 
generated by the Higgs mechanism), which under the assumption that it is sufficiently large has 
been chosen as an UV cutoff for the term (7). The reason for such a choice of the UV cutoff is 
based on the large distance behaviour of the bilocal correlator, described by Eqs. (15) and (16) 
below, according to which the mass m could be identified with the inverse correlation length of the 
vacuum in QCD, T~ x (see also discussion in Ref. 13). The latter quantity is large in the confining 
regime w.r.t. r -1 , where r is a typical length scale in the problem under study, i.e. a size of the 
Wilson loop in QCD, or a radius of the closed surface S in our approach. This fact served as a 
motivation for the curvature expansion of the string effective action in Ref. 6, which is actually 
an expansion in powers of (j^j . 

The propagator of the Kalb-Ramond field in Eq. (8) has the following form (the details of 
derivation of this propagator are presented in the Appendix to Ref. 13) 



D 



\u,pp 



X 



+ D 



(2) 

Xu,pp 



where 



(9) 



D 



(2) 

\i/,pp 



[x) = 



7/ 



Air 2 Nex 2 



— + - (K + K 2 

\x\ 2 



2 \x\ 



(10) 



' ^^Ap^/i^-i/ ^pyX\X p $p\X v Xp SypX^Xx 

Ki = Ki(m\x\),i = 0, 1, 2, 3, stand for the Macdonald functions, and one can show that the term 

J da X y{x) j dappix'jD^ix - x') 
s s 

could be rewritten as a boundary one, and therefore vanishes, since S is closed. In particular, 
as it could be anticipated from the beginning, when m — > 0, the contribution of the function (9) 
vanishes, whereas 
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J daxv(x) J da^x^Dfl^x - x') 



1 



2 (nNe)' 



dx„ & dx[ 



1 



^{x-x'y 



S S v ' dS dS 

Comparing Eqs. (8)-(10) with the result of the cumulant expansion of the weight factor (2), 
we obtain the following values of the coefficient functions, which parametrize the bilocal correlator 
of the dual field strength tensors 



D (mV) = 



171 

16V 2 



9- — - — me 



\x\ 



and 



D\ (m 2 x 2 ) = — 

V ' o7T 



9m 



2 X 2 



K\ m 



\x\ 2 

whose asymptotic behaviours at \x\ <C — and \x\ ^> — read 



K + K 2 



D 



9 m 2 



16tt 2 



x/ 



9 



47T 2 \x\ 



and 



D 



D 1 



9m 4 



16V2TT 
9 m 4 



m \x\ 



-m\x\ 



(m \x\ 



(11) 
(12) 

(13) 
(14) 

(15) 
(16) 



8V2n^ ( r 
respectively. 

We see that according to Eq. (8) all the terms of the cumulant expansion, higher than the 
quadratic one, vanish, so that the bilocal approximation to MVC is an exact result in the Londons' 
limit, or in another language, the ensemble of fields in this case is exactly Gaussian. 

The asymptotic behaviours (13)-(16) have already been discussed and compared with the 
lattice data 14 in Ref. 13. Namely, we see that Eqs. (15) and (16) reproduce correctly the long- 
range dynamics, i.e. the exponential fall-offs, and the result, that Di at large distances is smaller 
than D. However, as far as the short- distance dynamics is concerned, one can see that though Eq. 
(14) yields the correct ^r-behaviour of the function Di, Eq. (13) agrees with the results obtained 
within the MVC in QCD in Ref. 15 only in the lowest order of perturbation theory, where the 
function D is much smaller than the function D\. This result is a consequence of the Abelian 
character of the model under study. 



3. Realistic Values of Charge and V.E.V. of the Higgs Field and 
Consistency of the Londons' Limit 



The derivative (or curvature, or — -) expansion of the action 



6 



(iF) / da ^ X) I d ^^~ 4{X ~ X ' ? - ^ J d ^(x) J da, p (x')D^ p (x - x>) 
s s s s 

yields the following values of the string tension of the Nambu-Goto term and the inverse bare 
coupling constant of the rigidity term 

a = tit] 2 (9K (Ne) - 2) ~ ixif ^9 In - 2^J (17) 

and 

1 5tt 



«o ~ 8 (iVe) 2 



(18) 



respectively, where 7 ~ 1.781 stands for the Euler constant. These quantities are similar to the 
ones obtained for the low-energy expansion of the string effective action of the UCST in Ref. 10. 

We see that though the coupling constant (18) is always negative as it should be, the string 
tension (17) is positive only when the following inequality holds 

2 

Ne < — ~ 0.899, (19) 

7ea 

i.e. only for sufficiently small values of the Higgs field charge. At this point it looks natural to 
address the important question on how exact is the Londons' limit of AHM in modelling the real 
confining regime of QCD? In order to answer this question, one should identify the mass m with 
T" 1 , and substitute the value of the latter in the S'?7(3)-gluodynamics 14 , T~ x ~ lGeF, and the 
value of the string tension a ~ 0.2Gey 2 (see, for example, Ref. 5) into Eq. (17). This yields the 
following equation 

iVe = 1.12e- a22 ' a01 ( JVe ) 2 , 
whose solution according to MAPLE reads 

iVe ~ 0.892, (20) 

which luckly satisfies inequality (19). This means the consistency of the Londons' limit of AHM in 
modelling QCD in the confining regime. The corresponding square root of the Higgs v.e.v. reads 

r] ~ 1.121GeV. (21) 

Eqs. (20) and (21) provide us with the parameters of AHM in the Londons' limit, which according 
to the lattice data are the most adequate ones for the description of the confining phase of QCD. 

4. Equation of Motion of the Weight Factor of the String Represen- 
tation in the Loop Space 

In this short Section, we shall derive a loop equation for the weight factor (6) without the 
multiplier (7), which could be referred to the measure Vx^). We shall denote the weight factor 
(6) without this multiplier as (jF^ )■ 

Variation w.r.t. h uX (x) under the functional integral yields 
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^□iA + Wv + wJ -¥^h v ^Q[h[x($]]+Mitf vX *[h[x(0]]) =0, (22) 



where 



*[h[x(£)\\=exvUm J d<v(z(0)V MOl), 



and 

Eq. (22) when being reformulated in the loop space reads 



- ( iVe ) 2 ~x — ~~T~\ ~ l ^ 2 T uX ^(^l J [^((7)]) = 0, (23) 

where x M (cr),0 < <r < 1, is an element of the loop space, corresponding to an arbitrary closed 
contour C lieing on S [for example, x^(a) = x^ £ 2 = const.)], 

(7+0 





= J dr 



(7-0 

and 

+o 



&r M (r)' 



&V(x(<r)) ^ o ' 5^ ((7 + ±r) fa,, (<r - \t)' 



Eq. (23) coincides with the low-energy limit of the loop equation for the Wilson average in 4D 
compact QED, studied in Ref. 11, which is in the line with the discussion following after Eq. (7). 
The reader is referred to Ref. 11 for further investigations of Eq. (23). 

5. String Representation for the Partition Function of the Higgs 
Currents' Correlators 

Let us derive the bilocal correlator of the dual field strength tensors from the string represen- 
tation for the partition function of an arbitrary number of the Higgs currents. In the Londons' 
limit, such a partition function reads 



-\F^-^(d,e-NeA,) 2 + J,j, 



where = —Nerf{d ll 9 — NeA^) is the Higgs current in this limit. We use here the definition of 
Z [Jfj] normalized to unity, since in what follows we shall get the function D from the correlator 
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(jp( x i)j<r( x 2)) itself, rather than from the comparison of two exponents, as it has been done in 
Section 2. 

Performing the same duality transformation as the one used in Section 2, we get the following 
string representation for Z [ J M ] 



^-Jdx.J 2 (x) J da Xv (x) J da pp (x')D Xl/ ^ p (x-x') 

Z ^ = 7V—V—TZ — Vx ^ e 3 

V)£K„ „ Jx-x') J 



-n 2 J da Xv (x) J da pp (x')D Xl/ ^ p (x-x') 



exp^Nee X uaf3 J dxdy 



Ne 



d 2 



4 £ ^ 5 \dx a d Vl 
d 



D X u,np{x - y) )J(}{x)J 5 {y)+ 



+^T lip {y) I t^D Xu ^ p (x - y) \ J p {x) 



(24) 



Upon the two-fold variation of Eq. (24) w.r.t. J^, we arrive at the following expression for the 
correlator of two Higgs currents 



(NeY 



f dcr Xv (x) J d<T t _ lp (x')D Xv ^ p (x-x') 



fVx,(£)e s 

f da x „{x) J da pp {x')D Xv>pp (x-x') 



rfSp^xx - x 2 ) + 



1 d 2 



2dx la dx 2 , 



+n 2 j da 5Q (y) J da X9 (z)y-^—D XvM (x 1 - y) 







D„ 



> 2 - Z) 



(25) 



s s 

The first term in the figured brackets on the R.H.S. of Eq. (25) obviously vanishes for x\ ^ x 2 . 
It is also easy to show that the contribution of the term (10) to the R.H.S. of Eq. (25) vanishes 
too. As far as the part of the axionic propagator given by Eq. (9), is concerned, one can see that 
it has the order of magnitude of (Ne) 2 rf. Let us now assume that Nerf \S\ <C 1, where \S\ stands 
for the area of the surface S, so that the second term in square brackets on the R.H.S. of Eq. (25) 
is much smaller than the first term, and could be disregarded. This approximation corresponds 
to the case of a very small rj, i.e. to a very narrow potential. 

Taking all this into account and making use of the equation 5 , 



(j/3(xi)ja(x 2 )) 



d 2 



dx lu dx 2u 



<w - 



d 2 



dx la dx 2o 



D (m 2 (xi - x 2 ) 2 ) , 



which could be obtained from the classical equation of motion, we get the following (semiclassical) 
value of the function D 



2 2 

m x 



m 3 K\ (m \ x\) 

4"7T 2 \x\ 



, x^O. 



(26) 
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One can see that Eq. (26) differs in a factor | from the leading-order short-distance dynamics of 
Eq. (11), given by Eq. (13) (or from Eq. (11) itself at x ^ 0, where the last term on the R.H.S. 
without regularization vanishes). Besides the neglection of the second term in square brackets on 
the R.H.S. of Eq. (25), the difference between | and 1 describes the error of the semiclassical 
approximation. Notice also, that only the function D could be obtained from the correlator (25) 
in the semiclassical approximation due to the independence of the latter of the function D\ in this 
approximation. 

6. Conclusion 

Let us now stress some physical consequences of the obtained results and discuss open ques- 
tions. The first outcome of the Letter is the result of Section 2 that at large distances AHM in 
the Londons' limit coincides with the low-energy limit of the 4D compact QED. Secondly, we have 
demonstrated that the bilocal approximation to the MVC is exact in the Londons' limit, i.e. all 
the higher cumulants vanish. Thirdly, we have proved that the Londons' limit is consistent for 
the description of the confining phase of QCD according to the lattice data, and found the most 
adequate to QCD values of charge and the square root of the v.e.v. of the Higgs field, which are 
given by Eqs. (20) and (21) respectively. The most important result of the Letter is the two scalar 
functions, D and Di, given by Eqs. (11) and (12) respectively, which parametrize the correlator 
of the two dual field strength tensors and are known to be the key quantities in the MVC. 

Notice, that already in Ref. 13 it has been conjectured that the functions D and D 1 must 
follow from the propagator of the massive Kalb-Ramond field. In Ref. 13, this conjecture has 
been supported by comparison of the asymptotic behaviours of these functions at small and large 
distances obtained from such a propagator with the ones observed in the lattice measurements 14 . 
In this Letter, we have proved this conjecture and in addition obtained integration over metrics 
resulting from the integration over the singular part of the Higgs wave function. However, it is 
unclear what collective variable in QCD may serve for getting such an integration over metrics. 
As it has been already mentioned in the Introduction, the most physically appropriate candidate 
for this variable is the strong background fields, which provide confinement in the QCD vacuum. 
The solution of the problem of whether this is really the case or not will be the topic of a separate 
publication. 

It is also worth mentioning, that, as it has been discussed in Section 2, the short- distance 
asymptotics (13) of the function D agrees with the MVC only in the lowest order of perturbation 
theory. In order to get better correspondence with the MVC, one should investigate AHM in the 
vicinity of the Londons' limit, i.e. account for the corrections which come from the finiteness of 
the coupling constant A in Eq. (1), which presumably takes place for the "physical" AHM, i.e. the 
one, which describes QCD with the best accuracy 16 . In this way, the bilocal approximation will 
not be exact any more, and higher cumulants will start contributing into the Nambu-Goto term 
string tension and the rigidity term coupling constant. The effects arising due to these cumulants 
will be studied in the next publications. 

The final (technical) problem, which must be solved, is accounting for the Jacobian, which 
arises when one passes from the integration over the singular part of the Higgs wave function to 
the integration over metrics and has been omitted for simplicity in this Letter. This Jacobian has 
been calculated in Ref. 1 and, as it has been mentioned in the Introduction, yields the Polchinski- 
Strominger term 3 with the coupling constant ensuring exact cancellation of the conformal anomaly 
in D=4. 
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